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A fluid droplet located on a super-hydrophobic surface makes contact with the surface only at 
small isolated regions, and is mostly in contact with the surrounding air. As a result, a fluid in 
motion near such a surface experiences very low friction, and super-hydrophobic surfaces display 
strong drag-reduction in the laminar regime. Here we consider theoretically a super-hydrophobic 
surface composed of circular posts (so called fakir geometry) located on a planar rectangular lattice. 
Using a superposition of point forces with suitably spatially-dependent strength, we derive the 
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effective surface slip length for a planar shear flow on such a fakir surface as the solution to an 
infinite series of linear equations. In the asymptotic limit of small surface coverage by the posts, 
<f) s , the series can be interpreted as Riemann sums, and the slip length can be obtained analytically. 
For posts on a square lattice of periodicity L, our analytical results predict that in the low <j> s limit, 
£\J ■ the surface slip length, A, scales as 

-^ 3 / I" 3 , ,, /-. 

c: i~i6V^"^ ln(1+ ^ ) ' 

' O ■ which is in excellent quantitative agreement with previous numerical computations. 
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When a small droplet of liquid is deposited on a hydrophobic surface whose geometry is sufficiently rough, in certain 

<**> i situations the liquid will not fill the roughness grooves but will instead adopt a lower energy configuration, sometimes 

called a fakir state, where it remains on top of the surface topography (figure [1^). Exemplified by the Lotus leaf, 

surfaces for which this situation occurs are referred to as super-hydrophobic (or more generally, super-repellent), and 

their geometrical and physical description has been the subject of much recent work [lrl6| ■ 

One of the remarkable properties of super-hydrophobic surfaces is their low friction opposing the flow [7] . Since the 

I/"") ■ fluid next to such a surface makes contact with the solid at only a few isolated points and is mostly in contact with air, 

the shear stresses opposing fluid motion are small, and fluids in fakir states can move very easily. In these situations, 

the surface friction is usually quantified by a slip length, A, which is the fictitious distance below the surface where 

the no-slip condition would be valid on average J8l-ll0l|: no-slip corresponds to A = 0, while super-hydrophobic surfaces 

with A > have a lower effective friction than no-slip. 

From a practical standpoint, many different methods exist to design and produce super-hydrophobic surfaces [3|, |fjj, 

and they lead to surfaces with either random or controlled geometrical features. Surfaces with random topography 
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are known to show slip Ill4l3||. Planar surfaces with controlled geometry are of two kinds. The first kind are surfaces 



with one-dimensional features such as long grooves, usually aligned parallel or perpendicular to the flow direction, and 



for which a lot of experimental [14Hl9l | and modeling [18, 20 24] work has characterized their frictional properties. 

The second kind are surfaces with two-dimensional patterning, usually a series of vertical posts distributed on a 
regular lattice, as illustrated schematically in figure QJi. Only a handful of studies has focused on the friction of 



these types of surfaces, either experimentally [25( or theoretically 26|, |27|. In the limit of low solid fraction, these 



two-dimensional surfaces are expected however to display much lower friction (algebraic divergence of the slip length) 



than one-dimensional surfaces (logarithmic divergence) 26[, and it is thus important to be able to accurately predict 



their friction properties. This is the goal of the present paper. 
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FIG. 1: Fakir-like super-hydrophobic surfaces and setup for our calculation, (a): A droplet deposited on surface composed of 
tall hydrophobic posts can remain in a fakir super-hydrophobic state where it partially sits on the posts and partially on the 
air; (b): Definition sketches. A viscous fluid of shear viscosity fi sits in a fakir state on a two-dimensional array of circular 
posts. The posts have radius a, and are arranged on a rectangular lattice of periodicity D and L along the x and y directions 
respectively. A shear flow with shear rate 7 is imposed in the fluid along the x direction, and perpendicularly to the surface. 
In-between the posts, the fluid-air interface is assumed to remain flat and parallel to the super-hydrophobic surface. 



Previous theoretical work focusing on surfaces composed of posts on a lattice found that if <f) s < 1 denotes the 
areal density of the horizontal cross-section of the posts, the effective slip length of the superhydrophobic surface, 
non-dimensionalized by the typical distance, L, between the posts, behaves in the limit of low <f> s as 
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26|,|27|. In ([TJ, the positive coefficients, A and B, depend on the lattice and posts geometry but not 4> s , and can be 
fitted to numerical computations, while the square root scaling can be physically rationalized as follows. The shear 
stress acting on the posts scales as r ~ 4> s ^is where /z is the fluid viscosity and 7 S the typical shear rate around a 
post. If U is the typical fluid velocity between the posts, and a the typical post radius, we have j s ~ U/a, and since 

A ~ fJ-U/r we obtain A ~ a/4> s . As the solid fraction scales as (f> s ~ (a/L) 2 , we have a ~ L(f>J and hence we obtain 
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the square-root scaling of (UJ), namely A/L ~ l/<?v [2 

In this work, we consider the simplest possible geometry for flow over a fakir-like super-hydrophobic surface, namely 
vertical posts with circular cross section located on a regular rectangular lattice, and we provide two results. Firstly, 
we present an analytical method based on a linear superposition of flow singularities to accurately determine the 
friction (i.e. the effective surface slip length) for all surface coverage. Secondly, by asymptotically considering the case 
of low solid fraction <p s , we mathematically derive the scaling coefficients A and B governing (JT|), thereby predicting 
analytically the effective surface slip length. Our results are compared to previous computational work, and show 
excellent agreement. 



II. CALCULATION OF THE EFFECTIVE SURFACE SLIP LENGTH 

A. Problem outline 

Our geometry is illustrated in figureQJ). A viscous fluid of viscosity \i is located in a fakir state on a two-dimensional 
array of circular posts. The posts have radius a, and are arranged on a rectangular lattice of periodicity D and L 
along the x and y directions respectively, with la < min(£>, L). A shear flow with shear rate 7 is imposed in the far 
held by the prescribed velocity along the x direction, with shear perpendicularly to the surface along the z direction. 
The prescribed velocity is thus given by 

v ~ (jz + U)e x as z — > 00, (2) 

where e x , e y , e z denote unit vectors parallel to the axes, and thus U denotes the effective slip velocity at the fakir 
contact plane, z = 0. In-between the posts, the fluid-air interface is assumed to remain flat and parallel to the 
underlying super-hydrophobic surface with no shear acting on it. The flow satisfies the no-slip boundary condition on 
the top of the circular posts. 

Our aim is to accurately estimate the effective slip length, that is, the mean speed/mean shear ratio at the super- 
hydrophobic plane. As in our recent calculation for a moving grid |28| . the basic tool here is the Stokes flow generated 
by a similarly periodic array of aligned, parallel point forces directed in their plane. In addition, the force density 
function is assumed, for the leading approximation, to have an inverse square root rim singularity, as in the creeping 
flow past a disk. The resulting linear system has matrix elements that are interpreted as Riemann sums and the limit 
double integrals are easily evaluated. The dominant error arises from the missing term in the double sum and thus 
precise estimates can be obtained for the two coefficients in (p}. 

B. Calculation using a superposition of singularities 

The Reynolds number of the viscous incompressible flow is assumed to be sufficiently small for the velocity field v 
to satisfy the creeping flow (Stokes) equations [29J 

M V 2 v = Vp, V • v = 0, (3) 

where p is the dynamic pressure. The fluid motion can be represented as due to a distribution of tangentially directed 
Stokeslets over the circular regions, augmented by the uniform flow Ue x . The density functions must be both periodic 
in two dimensions and symmetric with respect to the sides of each rectangle. The field due to a two-dimensional 
rectangular array, periods D and L, of point forces of strength 4ir/iU\/DL directed parallel to the flow at infinity, is 
governed by 

V • v A = 0, (4) 
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The m\ — = ni2 term in © yields the anticipated shear at infinity and, as in [28], the flow generated by the 
oscillatory forcing is readily found by Fourier transform techniques [30]. Thus the flow governed by © is compactly 
expressed as 



VA = C/ 
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where, with a prime denoting that the mi = = 1712 term is missing, 
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with 5*i = V 2 S 2 - Only the velocities at the screen and at infinity are needed for the subsequent analysis. The solution 
(0 shows that 
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where, after substitution of flSJ) and ([3]), 



C(mi,m 2 ) = 1- -77^ 
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With suitably scaled force density functions in the circular regions, the total flow field is given by 

fn(a) 



v = Ue x + 
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in which the symmetry conditions restrict the density function to even cosines. The flow at infinity is determined by 
substitution of (flU)) into (fTU|) . which gives 



v- Ue x 
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a 2 VDL j -tv jo 
Comparison with ([5} then shows that the shear rate is given by 



Y^ fn(a) cos 2nj3 



adad/3, as z — > 00. 
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Thus only the mean force density contributes to the required slip length, A, given by 
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No-slip on the surface of the posts is achieved by enforcing v = at z = 0, r < a. The substitution of pT|) into 
([T3f gives the condition 
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where we have defined 

We next make two uses of the integral 
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cos2n/3exp[±27riMacos(/3 - ^)}df3 = 2(-l)"J 2n (27rMa) cos2m/', (19) 

where J p refers to the Bessel function of the first kind of order p, first for a direct evaluation in (TITl) which gives, 
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and then to write down the Fourier coefficients in (|2"0")) . Thus we obtain 

2 



<5fco = 



. f Y: M*)(-ir- k £ E ^^ m2 ) j 2w(27 rMa) 



a 2 VDL „,, 

n— mi m,2 



xJ 2fe (27rMr) cos 2ktp cos 2mp ada (0 < r < a, fc > 0), (21) 



where <5fc„ denotes the Kronecker delta. 
We next make two uses of the integral 
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first by introducing the conveniently scaled, inverse square root approximations for the force density functions, 



f n (a) ; (n > 0), (23) 



whose substitution in ([21]) yields 
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Then the immaterial r-dependence can be jettisoned by applying the operator, 

/• Q / r x2fc+i dr 

to (1241). in order to obtain 
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Equation (|25[) is a symmetric infinite system of linear equations for the coefficients {c n ;n > 0}. The substitution of 
into (jTHJ) then leads to the slip length as 

A=^- (26) 
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C. Asymptotic estimate of A in the low-<^) s limit 



An asymptotic estimate of the coefficients {c n ; n > 0} is found by noting that 
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is a Riemann sum, with increments (a/D,a/L) in (x,y). In terms of polar coordinates (r,ip), the integral has the 
separable form 
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where Neumann's symbol eo = 1, e n = 2 (n > 0). The corresponding solution of (|25|) is 
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whose substitution into (J16I) yields the dimensionless slip coefficient 
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where <f> s = tto 2 /DL is the fractional surface area covered by the circular regions. We note that this leading-order 
term is symmetric in (D,L), the respective periods along and across the imposed shear flow, corresponding therefore 
to an isotropic surface friction at this order. 

The leading error in the estimate pip is due to the absence of a term associated with mi = = m 2 in the Riemann 
sum. If the latter is regarded as having central function values, then the integrable singularity at the origin can be 
handled by exact integration over the rectangle centred at (0, 0), assisted by small argument estimates of the Bessel 
functions. Then (l29l) has the additional term 
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which is negligible unless n = k = 0. Evaluation for this case shows that (I5T1) has the more accurate form, 
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As a difference with the leading-order term, the next-order term in (|33[) is not symmetric in (D,L), and therefore 
rectangular lattices display anisotropic friction. Note also that this next-order term does not depend on the post 
radius a (and therefore <fi s ). For a square array (L = D), (|33j) simplifies to 
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FIG. 2: Results of our model, and comparison with previous computational and experimental work, in the case of a square 
lattice. Effective slip length, X/L, as a function of the solid surface fraction, <j> s (left), and as a function of l/y^> s (right). Solid 
line: Simple asymptotic low <j>g mode l. (1341) ; Empty squares: Computations by 27]; Filled squares: Computations by [26| ]; 
Filled circles: Experimental results by [251 ]. 



III. RESULTS AND COMPARISON WITH PREVIOUS WORK 



We now compare the results of our asymptotic model to previous work considering circular posts organized on a 
square array. The comparison is displayed in figure [21 where we show the effective slip length, X/L, as a function 
of the solid surface fraction, 4> s (left) and as a function of l/\f^ s (right). The solid line is our model, namely (|34l) . 
The symbols are data from previous work. Specifically, empty squares are computations by [27[, filled squares are 
computations by |26J, and filled circles are experimental results by |25| . 

The quantitative agreement between our model and previous numerical work is remarkable. The square root 



dependance on 
predict here A 
whereas 



evident in figure [5] (right) is indeed reproduce d by our model. For the coefficient A in ([T]) we 



r-f . 



3-y/7r/16 « 0.332. In their numerical simulations, 26] fit (JTJ to their results and obtain A — 0.325, 
27] fit to A — 0.34, which are both within about two percent of our prediction. 

For the second coefficient in ([T]), we showed above that a term independent of <p s is indeed the next-order term in 
the asymptotic expansion, and we obtained B — 31n(l + \/2)/2it ~ 0.421. [26j fit B = 0.44, which is less than five 
percent larger, while 27| obtain A = 0.468, a value about ten percent above our result. Overall, the integrated error 
between our simple model, (|34[) . and the numerics from [271 ] is about 1.8%, while the error between our model and 
the computations of [26] is about 3.9%. 

Our geometry implies that <j) s < 7r/4 but the asymptotic estimate (|34p loses validity at lower values, as seen on the the 
right-hand-side of figure [2j left), and as illustrated by its prediction of no-slip (A = 0) at <f> s = 7r 3 /[81n(l + v / 2)] 2 ~ 62%. 

Regarding the comparison with the experiments of [25|, we see that our model (and past computations) are able 
to predict the correct order of magnitude of the slip length, but do not capture the exact dependance on the surface 
friction, which appears to be weaker than inverse square root (see figure [5J right). This disagreement is likely due to 
the difference in flow field. In their paper, [25j measure the friction by using a cone-and-plate rheometer with rotation 
axis perpendicular to the surface place. As a result, the flow in their device has circular streamlines, and is not always 
aligned with the lattice periodicity. As a difference, in our work, as well as in [26j and 
planar shear flow aligned with the lattice period. 



27], the flow considered is a 



IV. CONCLUSION 



In this work we considered theoretically shear flow past fakir-like super-hydrophobic surfaces composed of circular 
posts located on a doubly-periodic rectangular lattice. Using a superposition of point forces with suitably spatially- 



8 

dependent strength, we obtained the total flow as an infinite series of linear equations. In the asymptotic limit 
of small surface coverage by the posts, the series can be interpreted as Riemann sums, which allowed us to derive 
analytically the effective surface slip length in the form of (p}. In the case of a square lattice, our results were found to 
be in excellent quantitative agreement with previous numerical computations. Future work will focus on embedding 
such low- friction fakir surface on a curved substrate (for example, a sphere), and attempting to predict its mobility 
coefficient in a viscous fluid. 



Acknowledgments 



We thank Chiu-On Ng and Christophe Ybert for providing us the numerical results from [27J and [261 re spectively, 
as well as Choongyeop Lee and Chang- Jin "CJ" Kim for providing us the experimental results from [25|. We also 
thank Lyderic Bocquct for useful discussions. This work was supported in part by the National Science Foundation 
(Grant CBET-0746285 to E.L.). 



[1] 
[2] 

[3; 
[4 

[5] 

[6] 
[7] 

[8] 

[9; 

[10 

[11 
[12 

[13 

[14; 

[15 
[16 

[17 
[18 
[19 

[20 

[21 
[22 
[23 
[24 
[25 
[26 
[27 
[28 
[29 
[30 



T. Onda, S. Shibuichi, N. Satoh, and K. Tsujii, Langmuir 12, 2125 (1996). 

J. Bico, C. Marzolin, and D. Quere, Europhys. Lett. 47, 220 (1999). 

L. Feng, S. H. Li, Y. S. Li, H. J. Li, L. J. Zhang, J. Zhai, Y. L. Song, B. Q. Liu, L. Jiang, and D. B. Zhu, Adv. Materials 

14, 1857 (2002). 

P.-G. de Gennes F. Brochard-Wyart and D. Quere, Capillarity and Wetting Phenomena: Drops, Bubbles, Pearls, Waves 

(Springer, New- York, 2004). 

D. Quere, Annu. Rev. Fluid Mech. 38, 71 (2008). 

P. Roach, N. Shirtcliffe, and M. Newton, Soft Matter 4, 224 (2008). 

J. P. Rothstein, Annu. Rev. Fluid Mech. 42, 89 (2010). 

C. Neto, D. R. Evans, E. Bonaccurso, H.-J. Butt, and V. S. J. Craig, Rep. Prog. Phys. 68, 2859 (2005). 

L. Bocquet and J.-L. Barrat, Soft Matter 3, 685 (2007). 

E. Lauga, M. P. Brenner, and H. A. Stone, in Handbook of Experimental Fluid Dynamics, edited by A. Y. C. Tropea and 
J. F. Foss (Springer, New- York, 2007). 

S. Gogte, P. Vorobieff, R. Truesdell, A. Mammoli, F. van Swol, P. Shah, and C. J. Brinker, Phys. Fluids 17, 051701 (2005). 

C. H. Choi and C. J. Kim, Phys. Rev. Lett. 96, 066001 (2006). 

P. Joseph, C. Cottin-Bizonne, J. M. Benoit, C. Ybert, C. Journet, P. Tabeling, and L. Bocquet, Phys. Rev. Lett. 97, 

156104 (2006). 

J. Ou, B. Perot, and J. P. Rothstein, Phys. Fluids 16, 4635 (2004). 

J. Ou and J. P. Rothstein, Phys. Fluids 17, 103606 (2005). 

R. Truesdell, A. Mammoli, P. Vorobieff, F. van Swol, and C. J. Brinker, Phys. Rev. Lett. 97, 044504 (2006). 

C. H. Choi, U. Ulmanella, J. Kim, C. M. Ho, and C. J. Kim, Phys. Fluids 18, 087105 (2006). 

D. Maynes, K. Jeffs, B. Woolford, and B. W. Webb, Phys. Fluids 19, 093603 (2007). 

P. C. Tsai, A. M. Peters, C. Pirat, M. Wessling, R. G. H. Lammertink, and D. Lohse, Phys. Fluids 21, 112002 (2009). 

E. Lauga and H. A. Stone, J. Fluid Mech. 489, 55 (2003). 

C. Cottin-Bizonne, C. Barentin, E. Charlaix, L. Bocquet, and J. L. Barrat, European Phys. J. E 15, 427 (2004). 

J. Davies, D. Maynes, B. W. Webb, and B. Woolford, Phys. Fluids 18, 087110 (2006). 

C.-O. Ng and C. Y. Wang, Phys. Fluids 21, 013602 (2009). 

C. J. Teo and B. C. Khoo, Microfluid. and Nanofluid. 7, 353 (2009). 

C. Lee, C. H. Choi, and C. J. Kim, Phys. Rev. Lett. 101, 064501 (2008). 

C. Ybert, C. Barentin, C. Cottin-Bizonne, P. Joseph, and L. Bocquet, Phys. Fluids 19, 123601 (2007). 

C.-O. Ng and C. Y. Wang, Microfluid. and Nanofluid. 8, 361 (2009). 

A. M. J. Davis and E. Lauga, Phys. Fluids 21, 113101 (2009). 

J. Happel and H. Brenner, Low Reynolds Number Hydrodynamics (Prentice Hall, Englewood Cliffs, NJ, 1965). 

H. Hasimoto, J. Fluid Mech. 5, 317 (1959). 



